Influence of carrier-carrier and carrier-phonon correlations on optical absorption and 

gain in quantum-dot systems 
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A microscopic theory is used to study the optical properties of semiconductor quantum dots. 
The dephasing of a coherent excitation and line-shifts of the interband transitions due to carrier- 
carrier Coulomb interaction and carrier-phonon interaction are determined from a quantum kinetic 
treatment of correlation processes. We investigate the density dependence of both mechanisms and 
clarify the importance of various dephasing channels involving the localized and delocalized states 
of the system. 
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I. INTRODUCTION 

In recent years, semiconductor quantum dots (QDs) 
have been studied extensively due to possible applica- 
tions in optoelectronic devices like LEDs, lasers, or am- 
plifiers 0,0. In the rapidly emerging field of quantum 
information technology, QDs have been successfully used 
to demonstrate the generation of single photons or cor- 
related photon pairs 0, 0, |1J] . Furthermore, the strong 
coupling regime for QD emitters in optical microcavi- 
ties has been demonstrated 0, • A common aspect in 
fundamental studies and practical applications of QDs is 
the critical role of dephasing processes. They determine 
the homogeneous linewidth of the QD resonances, limit 
the coherence properties of QD lasers and their ultrafast 
emission dynamics, and have a strong influence on coher- 
ent optical nonlinearities. Moreover, dephasing processes 
are intimately linked to lineshifts of the QD resonances. 

Optical studies of QDs have been recently focused on 
self- assembled systems which are typically grown in the 
Stranski-Krastanoff mode. The resulting QDs are ran- 
domly distributed on a two-dimensional wetting layer 
(WL). The interplay of carriers in the localized QD and 
delocalized WL states is mediated by the Coulomb inter- 
action of carriers in addition to carrier-phonon interac- 
tion. 

Theoretical studies of absorption and gain in QDs 
have been accomplished on the basis of multilevel optical 
Bloch equations @, or semiconductor Bloch equations 
(SBE) with screened exchange and Coulomb hole contri- 
butions 0- However, in these approaches the dephas- 
ing is merely a parameter entering the calculation via a 
T2-time. In Ref. ^(j the dephasing due to Coulomb in- 
teraction has been calculated using SBE with correlation 
contributions due to Auger-like WL-assisted capture pro- 
cesses (see below) where scattering integrals have been 
evaluated in terms of free-carrier energies. A more elab- 
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orate analysis of dephasing due to Coulomb interaction 
in quantum wells (QWs), which included non-Markovian 
scattering integrals based on renormalized ene rgie s, re- 
vealed quantitative modifications of the results . We 
show that for QD systems the situation is different due to 
the appearance of localized states with a discrete spec- 
trum. The calculation of scattering integrals in terms 
of free-carrier energies breaks down for processes due to 
Coulomb interaction which involve only localized states. 
It turns out that these are among the dominant processes 
for QDs containing more than one confined shell. For 
other processes involving localized states, their energy 
renormalization is also of enhanced importance. In the 
first part of this paper we present a non-Markovian treat- 
ment of dephasing due to Coulomb interaction where en- 
ergy renormalizations under high-density conditions typ- 
ical for QD lasers are self-consistently included and the 
importance of various scattering channels is analyzed. 

For the interaction of QD carriers with LA-phonons, 
extensive work has been devoted to the low tempera- 
ture regime, 0, 0, It is generally acknowledged 
that acoustic phonons dominate the dephasing at low 
temperatures. Nevertheless, several low temperature PL 
measurements [lti Il7| revealed LO-phonon replicas. In 
Ref. it was suggested that for temperatures above 
100K the dephasing due to carrier-LO-phonon interac- 
tion is of growing importance. The interaction of QD 
carriers with LO phonons is strongly influenced by hy- 
bridization effects 01 > which require the application of 
the polaron picture. A quantum-kinetic description of 
carrier-phonon scattering based on polarons has recently 
been used to explain ultrafast carrier capture and re- 
laxation processes in QDs 0|. In the second part of 
this paper, we present the corresponding treatment of 
dephasing processes. Then the combined influence of 
Coulomb interaction and carrier-phonon interaction on 
optical absorption and gain spectra is analyzed for quasi- 
equilibrium excitation conditions. 

In Refs. [2(j and [2l| PL spectroscopy measurements 
revealed homogeneous linewidths of several meV at room 
temperature. Comparable results have also been found 
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by four- wave mixing experiments [15j. Our microscopic 
calculation can reproduce these experimental findings 
even though a direct quantitative comparison, which 
is not the purpose of this paper, would require better 
knowledge of the specific QD parameters as well as pos- 
sible improvements of our QD model. 

The primary focus of this work is on QD optical spectra 
in the room-temperature high-density regime relevant for 
laser applications. For this situation, a quantum kinetic 
theory for carrier-carrier interaction and carrier scatter- 
ing with LO-phonons is developed. It is shown that non- 
Markovian effects and energy renormalizations play an 
essential role in QD systems. The theory is evaluated to 
compare the contributions of various scattering channels 
to dephasing and line-shifts as well as their excitation 
density dependence. Studies of the temperature depen- 
dence are restricted to T > 150K where the inclusion of 
LA-phonons is not necessary. 

Our investigations show that the efficiency of dephas- 
ing processes depends strongly on the carrier density in 
the system. For low carrier densities and room tempera- 
ture, the electron-LO-phonon interaction is the dominant 
mechanism, leading to the appearance of additional side- 
peaks in the optical spectra due to polaronic hybridiza- 
tion effects. Even though for high carrier densities the 
Coulomb interaction becomes the dominant mechanism, 
the electron-LO-phonon interaction remains important 
as it continues to contribute to the linewidth and line- 
shape of the QD resonances. 

For typical InGaAs QD parameters, a QD density 
of 10 10 cm~ 2 on the WL and carrier densities around 
10 10 cm~ 2 at room temperature, we obtain a homoge- 
neous linewidth of the ground state resonance with a full 
width at half maximum (FWHM) of 3.2meV in agree- 
ment with the experimental findings of Ref. [2C|. For 
lower carrier densities and decreasing temperature, we 
can reproduce the observed dephasing reduction while 
for elevated carrier densities and decreasing temperature, 
the reduction of the dephasing due to LO-phonons is 
compensated by increasing dephasing due to Coulomb 
scattering. 



II. THEORY OF OPTICAL ABSORPTION & 
GAIN CALCULATIONS 

The coherent optical excitation of semiconductors can 
be described in terms of interband transition amplitudes 
"0a = (e Q /i Q ) and population functions f a = (e^ea) and 
fa = (hah a ) for electrons and holes, respectively. For 
the annihilation and creation operators of electrons and 
holes, a combines the quantum numbers for cither WL 
or QD states. The dynamics of ip a (t) and fa' h {t) m 
response to the optical field E{t) is determined by the 
SBE which can be used to incorporate many-body 
interaction effects to describe dephasing, energy renor- 
malization and screening as demonstrated for QWs in 
Ref. [2i|. The optical absorption or gain coefficient can 



be defined via the response of the system to a weak opti- 
cal field E{oj) in terms of the linear optical susceptibility 
x(^) = P(i>o) / E{uj). Here P(uj) is the Fourier transform 
of the macroscopic polarization P(t) — ^ Q d* a ^ a with 
the interband dipole matrix element d a . 

A solution linear in the optical probe field E(t) im- 
plies that changes of the electron and hole population 
fa' h due to E(t) can be neglected. In the following 
we consider quasi-equilibrium distributions for electron 
and holes which are approximately realized in injection- 
current driven devices. Alternatively the carriers can be 
generated with an optical pre-pulse where a sufficient de- 
lay of the probe pulse ensures the the excited carriers are 
thermalized and that coherent effects of the pump pulse 
can be neglected. 

In this regime, the Fourier transform of the coherent 
interband transition amplitude ip a obeys the equation 



(*» - <' HF - # HF >«(M + [!-/«- fi] (M 
= s£ oul (M + sl hon (hu) . 

(i) 

The single particle energies entering Eq. Q in- 
clude already renormalizations due to Hartree-Fock (HF) 
Coulomb effects according to e^ ,HF = e a a + S°' HF with 
a = e,h. Here e a a are the free-carrier energies and 
S a HF = Y,fi( v at3f3a - V a p a p) ffi combines the Hartree 
and exchange Coulomb self-energy. In a spatially ho- 
mogeneous system like bulk semiconductors or for the 
in-plane motion for QWs, the Hartree terms cancel due 
to local charge neutrality. In the QD-WL system, the ab- 
sence of local charge neutrality leads to Hartree contribu- 
tions for the QD states. The influence of the WL states 
can be incorporated approximately within screened QD 
Hartree contributions as discussed in detail in Ref. |24| . 
The Coulomb exchange contributions to the interband 
Rabi energy fl^ F (huj) = d ■ E{ftw) + J^pVapapippito) 
give rise to the excitonic resonance of the WL and QD 
states. 

Without the terms S° oul (tkj) and S^ hon (huj) we would 
recover the well-known semiconductor Bloch equations 
|22| formulated in an arbitrary eigenfunction basis. 
S2° ul (fiuj) and S Fhon (fkjj) are used to describe correlation 
contributions which lead to dephasing of the coherent 
polarization and to the corresponding additional energy 
renormalizations due to Coulomb interaction and carrier- 
phonon interaction, respectively. These terms will be dis- 
cussed in the following sections. A general feature of the 
following description is that the dephasing contributions 
for both the Coulomb interaction and the carrier-phonon 
interaction can be separated into a part which is diagonal 
in the state index a for the transition amplitude (diag- 
onal dephasing r DD ) and a corresponding off-diagonal 
part (r OD ) which mixes the coherent interband transi- 
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tion amplitudes for various states according to 

s a (hu) = -r° D (ftw)v> a (M 

+ E r °£(MV> Q1 (M- (2) 

Oil 

The imaginary part of r DD describes a dephasing (damp- 
ing) of the coherent polarization which can be expressed 
as a state and frequency depcndend T 2 time. The con- 
tribution of r DD is however partly compensated by the 
off-diagonal contribution r OD . The real parts of r DD 
and r OD give rise to additional renormalizations of the 
energies on the L.H.S. of Eq. 0). 



III. COULOMB INTERACTION 



A. Theory 



The correlation contributions due to Coulomb interac- 
tion, S~ oui {tkj), are evaluated in the second-order Born 
approximation. To close the set of equation for the in- 
terband transition amplitude ip we additionall y apply 
the generalized Kadanoff-Baym ansatz (GKBA) |25Ll26| . 
Within this approach the correlation contributions due 
to Coulomb interaction read 



rT(M = i E E 

a,b—e,h 0.10.20.3 

x{ W aa . 2Cl3ai [2W* a2asai -W* a 

+ 2W cea3Q , 3ai W aa2 a 3 oii 



, ia „] 9 {hw - s b a - e a ai + (e a a2 )* - e« 3 ) [ (l - / Q a 2 ) f a a J^ + (/ - 1 - /) 
g (hw 



it 



^-^ 2 + (4 3 ) /a 2 (l-/a 3 ) 



and 



= * E E 

a,b=e,h CK2C*3 
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1 - 



(3) 



/)]} , 



(4) 



: { W aa2asai [2W* a2a3ai - W: a2aia3 ] g (tkJ ~ i a a - e b ai - e a a2 + (e a a3 )*) [ (l - /£,) f a a J a a + (/ -> 1 - /) 
f 2^ aQ2Q3Ql TT ( : Q2Q3Qi .9 {hw - e a a - e b ai + (s b J* - e b a3 ) [f b 3 (l - f b J /« + (/ -> 1 - /)] } , 



with the screened Coulomb interaction matrix elements 
W aa2a3ai discussed in Appendix ^ Note that e° = 
e a — appearing in the function 5(A) = are ef- 
fective complex single-particle energies which combine 
renormalized energies as well as the corresponding 
quasi-particle damping 7*. This is elaborated in more 
detail below. 

As for the particle scattering [2]} we have direct 
and exchange contributions which are proportional to 
2\W aa2 W* a2aia3 , respectively. The 

population factors describe the availability of initial and 
final states. For the above discussed excitation condi- 
tions, the population factors are time independent and 
we have restricted ourselves to the contributions linear 
in the transition amplitude ip a . 

The frequency dependence of r DD and r OD in Eqs. © 
and J2J reflect the non-Markovian treatment of the de- 
phasing processes. Indeed, in this case, the product of 
the type Tip of Eq. (0) amounts, in the time domain, to 
a convolution integral describing memory effects. The 
Markovian limit is obtained by pulling out from this in- 
tegral the slowly varying component of ip a (t) . Then the 



fast component e^ £ <* +e <»)* (a / b) fixes the frequency 
value at ui = ui a = + £„) for r DD and uj — u ai for 
r OD , so that in this limit both dephasing contributions 
become frequency independent (see Appendix [U] for de- 
tails). If one combines the Markov approximation with 
the use of free-carrier energies in the scattering contribu- 
tions, the diagonal dephasing is given by the sum of in- 
and out-scattering rates, Im(r° D ) = + S° ut . Here 
gm.out are .j-kg ra t es j n th e Markovian kinetic equation 
for the carrier population, as defined e.g. in Ref. [27| . 
Results of these equations for optical spectra have been 
studied in detail for QW systems. If one restricts the 
analysis to diagonal dephasing contributions r DD and ne- 
glects off-diagonal dephasing terms r OD , dam ping of the 
excitonic resonances is grossly overestimated 23]. Non- 
Markovian calculations further reduce the interaction- 
induced broadening and line-shift of the excitonic res- 
onances in QW systems [Tlj. 

For the following discussion we would like to point out 
that one has to compare on one side the Markov approxi- 
mation with a non-Markovian treatment and on the other 
side the use of free-carrier energies versus renormalized 
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energies in the scattering integrals. The non-Markovian 
treatment is more crucial in QD systems compared to 
QWs, since the discrete part of the spectrum with large 
energy separation emphasizes the frequency-dependence 
of the dephasing rates [T(j • 

For the use of free-carrier energies in the scatter- 
ing contributions, the limit 7^ — > 0, 7^ > leads to 
g(A) = n8(A)+iP-g where P denotes the principal value 
integral. This approximation implies serious difficulties 
in a QD system when processes are taken into account 
which involve only discrete states. In these cases the 
^-functions are not integrated out and thus the results 
are not well defined. Note that this is not a artifact of 
the Markov approximation but also applies to the non- 
Markovian calculation. Even if one would introduce a fi- 
nite broadening of the S- function, e.g., due to interaction 
of carriers with acoustic phonons, the non-Markovian cal- 
culation still yields unphysical results. This is because 
the spectrum predicted by Eq. has peaks at the ener- 
gies of the L.H.S. renormalized by the correlation contri- 
butions on the R.H.S. and hence, hu> samples renormal- 
ized interband transitions. If they are mixed with free- 
carrier energies in the correlation contributions Eqs. 
and via g(fruj — A), then artificially also free carrier 
transitions and hybridization effects between all these en- 
ergies can appear in the optical spectra, which are absent 
when self-consistently renormalized energies are used. 

Another important point, which complicates the treat- 
ment of the coupled QD-WL system, is the existence 
of many scattering channels which are often partly ne- 
glected (for a discussion of the scattering channels in the 
framework of carrier kinetics see Ref. [23). This can 
strongly influence the results as we will show in Section 

mm 

Self- consistent single-particle energy renormalizations 

The renormalization of the single-particle energies, 
which enter in the dephasing rates, originates from the 



same interaction processes that determine the dephasing 
itself. Technically, they can be traced back to the same 
many-body self-energy. The renormalized energies can 
be obtained from the poles of the retarded Green's func- 
tion (GF) in Fourier space. In the absence of interaction, 
the retarded GF is given by 



with the single-particle energy e a a and 6 — > 0, 6 > 0. At 
elevated carries densities, it is a reasonable approxima- 
tion, that the main effect of the Coulomb interaction is 
a shift of the single-particle energy and the addition of a 
quasi-particle broadening such that the single-pole struc- 
ture of the retarded GF remains valid. This picture cor- 
responds to the Landau theory of a Fermi liquid and leads 
to the retarded GF of the interacting system 



'f a (M = — - - = — . (6) 

' ' to-e% kuj~e a a + i^ y ' 



With this ansatz, we define a complex effective single- 
particle energy e° = e a a — ij® which consists of a 
renormalized energy e a a and the corresponding quasi- 
particle damping 7™. In the pole approximation, the 
self-consistency requirement leads to 



e a = e a + S Q ' + S Q (ej (7) 



with the retarded self-energy 



sr ct (M = -i E E 

b—e,h ct\a.20t.3 

x{ W* a2Ct3ai [2W* a2a3ai - W* a2aia3 ] 9 (fcj - e a ai + (i a a2 Y - e a a3 ) [ (l - /£) /£/« +(/-!-/)' 



(8) 



+ 2W QQ2Q;3Q1 W aa2a30tl 



9 fa - e« ai - 4 + (4J*) [fl (1 - fi 3 ) r ai + 



The close connection between the diagonal dephasing and 
the retarded self-energy £^ ret can be expressed as 



B. Results 



E 

,,b—e,h 



(to - i» a ) 



(9) Throughout this paper, results are presented for a car- 

rier and lattice temperature of 300K and a density of 
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QDs on the WL of 71qd = 10 10 cm -2 . Further details of 
the QD model are given in Appendix 1X1 

In Fig. Q] we show optical absorption spectra for the 
combined QD-WL system and various carrier densities. 
Excitation-induced dephasing and energy renormaliza- 
tions due to Coulomb interaction have been included ac- 
cording to Eqs. ©,0 and Identifiable are the 
excitonic resonance of the WL at around -15meV as well 
as the p-shell and s-shell resonances at about -90meV and 
-150meV, respectively. We observe a strong damping of 
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FIG. 1: Imaginary part of the optical susceptibility for the 
combined QD-WL system including interaction-induced de- 
phasing and line shifts due to Coulomb interaction for various 
total carrier densities. The inset shows a scale up of the QD 
resonances. 



these resonances with increasing carrier density which is 
accompanied by a pronounced red shift of the QD lines. 
The transition from absorption to gain takes place at a 
density ~lx 10 11 cm~ 2 (~5x 10 11 cm~ 2 ) for the s-shell 
(p-shell). An important point for practical applications 
is the saturation of the optical gain for the s-shell ac- 
companied by an increasing red-shift. This is due to a 
combination of state filling and saturation of dephasing 
at the s-shell resonance. In contrast the gain at the p- 
shell resonance increases further and shows no saturation 
for the densities investigated here. 



Importance of different dephasing processes and frequently 
used approximations 

The Coulomb matrix elements allow to distinguish be- 
tween different dephasing processes in complete anal- 
ogy to the carrier scattering If all four indices of 



processes as WL relaxation. They describe dephasing 
and energy shifts of the WL states in the optical spectra. 
Apart from the fact that we have to include OPW cor- 
rections to the interaction matrix elements for a proper 
description of the coupled QD-WL system [271 , this re- 
sembles the case of a QW system. Since the processes 
involving QD states (except intra-dot processes, see be- 
low) are discussed in great detail in Ref. we will only 
give a brief survey here. 

If one of the four indices of W aa2a3ai corresponds to 
a QD state, we consider the dephasing mechanism as 
WL assisted carrier-capture, because in the correspond- 
ing scattering process an electron or a hole is captured 
from the WL to the QD jS^. Likewise we refer to pro- 
cesses with two QD state indices as WL assisted relax- 
ation and to processes with three QD state indices as 
dot-assisted processes. Intra-dot processes are described 
when all four indices of W aa2a3ai correspond to QD 
states. These scattering events are not important in the 
carrier dynamics, because we only consider two confined 
shells for electrons and holes and therefore such scatter- 
ing processes cannot redistribute carriers. However, they 
clearly provide additional dephasing of the coherent po- 
larization. Note that the intra-dot processes are not the 
only events which produce dephasing without redistribu- 
tion of carriers. Such processes, in which two carriers 
switch their states, also appear in the WL assisted relax- 
ation and in the WL relaxation. 
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FIG. 2: Imaginary part of the optical susceptibility for the 
s-shell resonances at a carrier density of 2 x 10 11 cm -2 . The 
result including all considered Coulomb scattering processes 
(thick solid line) is compared to calculations where only cer- 
tain classes of processes are evaluated. 



Figure |21 shows that WL assisted relaxation (dash- 
dotted line) and intra-dot scattering (thin solid line), 



6 



which include the processes leading to dephasing with- 
out redistribution of carriers, are the most important 
contributions, while WL assisted capture (dotted line) 
and dot-assisted (dashed line) processes cause a rather 
small dephasing. However, since this picture can change 
with carrier density, it is important to test the influence 
of all dephasing channels. 

In Fig. |3 we compare the optical susceptibility from 
the full calculation (solid line) with a result where only 
the diagonal dephasing contributions were taken into ac- 
count (dotted lines). For the excitonic resonance of the 



damped out in this approximation. 

The influence of the exchange term in the dephasing 
rates, which is often disregarded, is also investigated in 
Fig. [3J Neglecting the exchange scattering in the dephas- 
ing contributions clearly overestimates the homogeneous 
linewidth by about 30 %, thus pointing out that the ex- 
change terms are almost as important as the direct term 
and should be included in the calculation. 

Renormalization of single-particle states 



n SYS = 2xlO U cm 2 



r DD (w) + r° D (ro), with exchange scattering 
r DD (w) , with exchange scattering 

r DE 

(o)) + r (to), without exchange scattering 
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FIG. 3: Imaginary part of the optical susceptibility for a 
carrier density of 2 x 10 11 cm -2 with full dephasing contri- 
butions (solid lines), with diagonal dephasing contributions 
(dotted lines), and with diagonal and off-diagonal dephasing 
but without the exchange scattering in both terms (dashed 
lines). 



WL, the absence of the off-diagonal dephasing leads to 
a overestimation of the linewidth by roughly a factor of 
two under the considered high-density conditions. Using 
only the diagonal dephasing contributions turns out to 
be a reasonable approximation for the lowest QD reso- 
nances while for the excited QD transition the lineshape 
is not fully reproduced. Off-diagonal dephasing contri- 
butions are less important for the QD resonances due 
to the rather large spectral separation between the QD 
states and the WL states and for the same reason off- 
diagonal contributions are stronger for the p-shcll than 
for the s-shell. 

Note that the foregoing discussion applies to the non- 
Mar kovian treatment. If one neglects the off-diagonal 
dephasing and uses the Markov approximation, again 
the dephasing is grossly overestimated when all relevant 
scattering processes are included. For the discussed QD 
model, transitions due to localized states are completely 



As discussed above, the single particle renormaliza- 
tions are of critical importance for the proper determina- 
tion of the dephasing rates. In Fig. 0Jt the single particle 
renormalizations due to Coulomb interaction are shown. 
The Fock and correlation contributions lead to a decrease 
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FIG. 4: Renormalized single-particle energies of the QD 
states for electrons and holes as a function of the total carrier 
density in the system (a) and corresponding single-particle 
broadening(b). 



of the single-particle energies with increasing carrier den- 
sity. The Hartree terms increase the energies for electrons 
and decrease the energies for holes (for a detailed discus- 
sion cf. (24|). This leads to slightly increasing (decreas- 
ing) single-particle energies for the electrons (holes) from 
low to intermediate carrier densities. For higher densi- 
ties, both electron and hole energies strongly decrease. 

With increasing carrier densities the single-particle 
broadening (Fig. 2t>) strongly grows up to intermediate 
carrier densities due to the larger availability of scatter- 
ing partners. For higher carrier density, the strong single- 
particle energy shifts lead to a reduction of the scattering 
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efficiency. Additionally at very high carrier densities the 
intra-dot processes are suppressed due to Pauli-blocking. 
These effects result in a strong decrease of the single- 
particle broadening. 

IV. CARRIER-PHONON INTERACTION 
A. Theory 

In this section, we compute correlation contributions 
in the optical spectra due to interaction of carriers with 
LO-phonons. In contrast to the frequently applied time- 



r° D (M 



and 

C D (M 



where the prefactor M^ = 47ra^=(/^Lo) 3 ^ 2 includes 
the polar coupling strength a, the reduced mass m and 
the phonon energy huj^o- In this equations we have in- 
troduced the abbreviation 

G;< b » = Jdr G*' ret (r) G h a ^(r) (12) 

which represents the Fourier transform of a product 
of two retarded polaronic GFs. Please note that the 
structure of this quantity can be traced back to the 
change from the c, v picture to the e, h picture. Through 
these functions polaronic renormalization effects such as 
phonon replicas and hybridization between the localized 
states [l9( are included in Eqs. IjlOfl and Ijlljl. Details for 
the calculation of the retarded polaronic GF are give in 
Appendix [B] 

B. Results 

The high-density spectra with dephasing due to 
Coulomb interaction exhibit only three resonances, 
namely s-shell, p-shell and the excitonic resonance of 
the WL. Absorption spectra with correlation contribu- 
tions due to interaction of carriers with LO phonons 



dependent perturbation theory, it has been pointed out 
in Ref. |l|| that the interaction involving localized QD 
carriers requires a description in the polaron picture. In 
Ref. |l9j . quasi-particle renormalizations described by a 
polaronic retarded GF have been used with in a quantum- 
kinetic theory to evaluate scattering processes and pop- 
ulations changes. We extent this treatment to the polar- 
ization dynamics to analyze the corresponding dephasing 
and interband-energy renormalizations. 

Within the random-phase approximation (RPA) and 
the GKBA, the correlation contributions due to LO- 
phonons are 



are shown in Fig. [S] Polaronic renormalizations of the 
single-particle states lead to a more complicated reso- 
nance structure for the interband transitions displayed in 
the inset of Fig. [S] One can identify phonon replicas and 
results of the hybridization of the single-particle states. 
For example, the s-shell resonances has a shoulder on 
the lower energetic side due to hybridization of the cor- 
responding electron state. Energetically above the res- 
onances of the s-shell and the p-shell several peaks due 
to phonon replicas and their hybridizations can be ob- 
served. The non-Lorentzian character of the lineshapes 
is even more pronounced as for the spectra with Coulomb 
interaction. On the other hand, the damping of the res- 
onances at higher densities remains weak. 



V. DEPHASING DUE TO CARRIER-PHONON 
AND COULOMB INTERACTION 

In this section, we investigate the combined influence 
of the carrier-carrier Coulomb scattering and the inter- 
action of carrier with LO phonons. It turns out that 
polaronic resonances in the single-particle spectral func- 
tion are strongly damped out due to Coulomb scattering 
of carriers even at low carrier densities, so that a pole 
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Af 2 f 
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FIG. 5: Imaginary part of the optical susceptibility for the 
combined QD-WL system including interaction-induced de- 
phasing and line shifts due to carrier-phonon interaction for 
various total carrier densities. The inset is a logarithmic plot 
of the QD resonances. 



approximation for the single-particle properties of the in- 
teracting system is reasonable. On the other hand, in- 
terband energy renormalizations and dephasing contribu- 
tions due to the interaction of carriers with LO phonons 
continue to be important. Since our goal is a calculation 
of optical spectra for the QD-WL system under the in- 
fluence of correlation effects at elevated carrier densities, 
we use for the description of single-particle properties a 
retarded GF obeying Eq. I|B1|) with the free carrier en- 
ergy in the L.H.S. replaced by e a a — e a a + A° — £7®, which 
is determined by Eqs. J7J) and (JSJ. In other words the po- 
laron is obtained by dressing with the phonon interaction 
not the free particles but the quasi-particles obtained by 
Coulomb renormalization. Rewriting the retarded GF as 



Gr(r) = 0S(r) e-n 



-UK 



(13) 



we separate in Ga( T ) the phonon renormalization effects 
which are, of course, influenced by the presence of the 
Coulomb interaction. (The equation obeyed by Q still 
contains A and 7.) The finite life-time of these quasi- 
particles produces in general sufficient damping to reduce 
the polaronic GF to a single-pole structure. This pole 
is used instead of e a a in Eq. Q). In this this way the 
Coulombian and the polaronic problems become coupled 
and have to be solved self-consistently. The iterative 
solution to this problem converges rapidly. Results of 
the spectral function of QD and WL states are shown in 
Fig. El Even for low carrier densities, polaronic structures 
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are strongly broadened as a result of the dominant role 
of damping due to Coulomb scattering. (For low carrier 

la 

densities the exponential decay e~~ T due to Coulomb 
interaction, which is superimposed to the polaronic func- 
tion Q% (r) in Eq. (|13l) , might somewhat overestimate the 
damping of the polaronic resonances. Nevertheless at 
higher carrier densities we expect a strong broadening of 
the polaron satellites.) 

In Fig. [7| calculated absorption spectra are shown, 
which include correlations due to carrier-carrier scat- 
tering Eqs. 10-0, and interaction with LO-phonons, 
Eqs. (|lL)fl - l|ll|) . both evaluated with self-consistently 
renormalized single-particle energies. As in the result for 
dephasing due to Coulomb interaction, Fig. ^ we see the 
bleaching and red-shift of the resonances due to many- 
body interactions and also the saturation of the s-shell 
gain. Although we observe that the Coulomb interaction 
is the clearly dominant dephasing mechanism for high 
carrier densities, we also infer from a comparison of Fig.[7| 
with Fig. that even in the gain regime the electron- 
phonon interaction gives rise to a clear increase in the 
dephasing. Nevertheless polaronic features are absent in 
the spectra, since the complicated multi-peak structure 
of the spectral function is completely damped out due to 
Coulomb effects. 

For intermediate carrier densities around 5 x 10 10 cm~ 2 
both types of interaction processes are equally important. 
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FIG. 7: Imaginary part of the optical susceptibility for 
the combined QD-WL system including interaction-induced 
dephasing and line shifts due to Coulomb interaction and 
carrier-phonon interaction for various total carrier densities. 
The inset shows a scale up of the QD resonances. 



Comparing the results in Fig. ISlone can conclude for our 
situation that taking only the Coulomb dephasing mech- 
anism into account underestimates the dephasing of the 
ground state transition by roughly a factor of two, while 
the damping of the WL is even dominated by carrier- 
phonon interaction. For higher carrier densities, however, 
this picture changes as can be seen in Fig.^l Using a car- 
rier density of 2 x 10 11 cm" 2 , the Coulomb interaction is 
clearly the dominant mechanism for the QD resonances. 
For the excitonic resonance of the WL, the two mech- 
anisms are equally important even at this rather high 
carrier density where we are already in the gain regime 
for the s-shell transition. 



Comparison with Experiments 

In most experiments, the emission from an ensemble of 
QDs is studied such that the additional inhomogeneous 
broadening of the QD resonances due to size and compo- 
sition fluctuations contributes. However, recent experi- 
ments |2jJ, E^l have been performed at room temperature 
on single QDs. In Fig. 1(a) of Ref. [2l| optical spectra are 
displayed for different excitation intensities. Although 
a direct comparison of the excitation dependence is not 
possible, due to uncertainties in the carrier densities gen- 
erated in those experiments, the general features of the 
spectra are identical. A clear distinction between the 
QD transitions and the excitonic resonance of the WL 
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is seen. The QD resonances show a density dependent 
bleaching and a pronounced red shift due to many-body 
correlations while the spectral position of the WL is al- 
most unchanged and the WL resonance is only bleached 
out. This means that the red-shift of the QD resonances 
cannot be attributed to band-gap shrinkage effects of the 
WL. 
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The observed homogeneous linewidth in this experi- 
ment is 8-13meV depending on the excitation intensity. 
This is larger than our findings, but one has to take into 
account that the QD, investigated in Ref. |2^|, seems to 
have three confined electronic shells which results in more 
dephasing channels. A better comparison is possible with 
the results of Ref. ;20j because two types of QDs are in- 
vestigated, where one type resembles our model in the re- 
spect that it is supposed to have two confined electronic 
shells. Regarding the s-shell, we can infer from Fig. 3 
of Ref. [2(| that the observed homogeneous linewidth is 
slightly above 3meV at room temperature which is simi- 
lar to our result of 3.2meV. The calculated value is prac- 
tically unchanged for a carrier density range from 5 x 10 8 
to 10 10 cm- 2 . 

For a better comparison with experiments, we have 
also calculated the temperature dependence of the ho- 
mogeneous linewidth in a temperature range in which 
the influence of LA-phonons remains small. As shown in 
Fig. 1101 we obtain for low carrier densities (where carrier- 
carrier scattering is sufficiently weak) the expected reduc- 
tion of dephasing with decreasing temperature as seen in 
Ref. [20]. (A quantitative comparison would require a 
better knowledge of QD parameters and phonon modes 
in the QD-WL system.) For the case of a higher excita- 
tion density of 10 10 cm~ 2 , the line-broadening is almost 
temperature independent due to a balancing of different 
scattering channels. The dephasing due to LO-phonons 
decreases with temperature as in the low-density case. 
However, the Coulomb interaction is no longer negligible 
for this carrier density. With decreasing temperature, the 
WL states are less populated while the occupation of QD 
states increases. This enhances intra-dot relaxation pro- 
cesses and provides stronger dephasing due to Coulomb 
interaction. 
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VI. CONCLUSION 

We have calculated the excitation induced dephasing 
and lineshifts for a QD-WL system on a microscopic 
basis. Both Coulomb and LO-phonon contributions to 
the homogeneous linewidths are found to be equally im- 
portant for elevated carrier densities. The role of self- 
consistent single-particle energy renormalizations in the 
scattering integrals is emphasized. For the Coulomb in- 
teraction the relative importance of various scattering 
channels has been analyzed. 
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APPENDIX A: QUANTUM DOT MODEL 
SYSTEM AND COULOMB MATRIX ELEMENTS 

We focus on lens-shaped InGaAs QDs which are lo- 
cated on top of a two dimensional WL. Such systems typ- 
ically consist of localized QD states energetically below 
a quasi-continuum of delocalized WL states for electrons 
and holes. In the present paper we assume the single 
particle spectrum as depicted in Fig. 1 of Ref. j23| • For 
the in plane motion we assume a parabolic confinement 
potential, leading to harmonic oscillator like single parti- 
cle wavefunctions. To account for the finite depth of the 
confinement potential, we take into account two confined 
shells for electrons and holes, which we refer to as s-shell 
for the ground state and p-shell for the double degener- 
ated excited state, due to their angular momentum prop- 
erties. We choose a level spacing of 40meV for electrons 
and 15meV for holes so that the free carrier transitions 
appear at -55meV relative to the WL bandedge for the 
p-shell and at -HOmeV for the s-shell. The WL states are 
modeled by orthogonalized plane waves (OPW) j2?|]. For 
the motion in growth direction (z-direction) we assume 
a infinite height potential well of 4 nm thickness. 

Under this assumptions we can construct the Coulomb 
matrix elements 



V, 



0.(3^8 



x / d 2 Q(p* a (Q)ifs(g)e 



iq Q 



x / <*Y (pp(g')^(e') e 



tq.g 



(Al) 



consisting of overlap integrals between single-particle 
wavefunctions and the Coulomb potential Vq as given 
in Ref. [27| • For the screening we use a generalization of 
the static Lindhard formula which is also explained in de- 
tail in Ref. [2?! ■ This procedure leads to the replacement 
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Vq — > Wq in Eq. (|A1|) for the matrix elements W a pjs- 
Because we are working in the envelope approximation 
with equal envelope wavefunctions, we do not need to 
consider band indices for the Coulomb matrix elements. 
All other parameters are chosen as in Ref. [27j. 



where we consider monochromatic LO-phonons with the 
frequency u>lo- The corresponding population of the 
phonon bath is given by a Bose-Einstein distribution 
riLo = l/(e riWLO / fcT — 1) and Vp a p a is the unscreened 
Coulomb matrix element (cf. section lA"|l . 



APPENDIX B: RETARDED POLARONIC 
GREEN'S FUNCTION 



The retarded polaronic GFs obeys the equation of mo- 
tion 



GT\r) = S(r) 



Jdr' E^ rct (r-r') G^ rct (r' 



(Bl) 



The corresponding retarded self-energy is given in RPA 
by 

S a / Ct (r) = ihJ2 G^ ct (r) D<J-r). (B2) 





APPENDIX C: CONNECTION BETWEEN 
MEMORY EFFECTS AND FREQUENCY 
DEPENDENCE IN THE SCATTERING 
INTEGRALS 



In this section we is show how the non-Markovian 
(Markovian) scattering integrals lead to frequency de- 
pendent (independent) dephasing terms. The time do- 
main formulation of the SBE with correlation contribu- 
tions due to Coulomb interaction is given by 



We assume for the calculations restricted to electron- 
phonon interaction that the polaronic retarded GF is not 
influenced by population effects. This has been verified 
for a bulk system in Ref. |2^| . 

Assuming that the phonon system is in thermal equi- 
librium, the phonon propagator (combined with the in- 
teraction matrix elements) is given by 



m d<(t) 



Mlo 
e 2 /e 



iVi 



lo e 



+ (i + tvlo) e 



(B3) 



[in- - eT" - )<Mt) + [i - r a - /£] o„ (t) 



= — ih s. 



Coul 



(*) 



(CI) 



In the following we explicitly consider the direct e-e and 
h-h interaction contributions of the diagonal dephasing, 



s a (t) = re J *> £ ]T 2W aa2Ct3ai (t)w* a2a3ai (t') 



x G^' ct (i, t') [G^\t, t')} * G*f *(*, t')G^ ct (t, t>) Mt')F a2a3ai (f) 

I 



(C2) 



Other terms can be treated in complete analogy. Here 
we have defined 



which corresponds to Eq. ©, we obtain 



+(/-»!-/)• (C3) 

We assume static screening where the Coulomb matrix 
elements depend only parametrically on time via the pop- 
ulation functions. For the discussed excitation conditions 
both / and W are time independent. With the ansatz 

G a a ict (t,f) = -ie(t_t') e -*n(*-f) , (C4) 



S a (t) — ^ 2 $Z 2W aa2a3ai W aa2a3ai F a2a:iai 



a,b—e,h aia 2 as 
b=£a 
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With the Fourier transform ip a (t) = J ^ e ~ iujt ip a (u) 
and using the integral relation 

t 

J dt'e* A(t -^ = J ImA > (C5) 

— oo 

we arrive at Eq. (JJJ and the terms of Eq. which are 
discussed here. 

For the Markov approximation one separates from the 
transition amplitude a phase factor, which rapidly oscil- 
lates with the (renormalized) interband transition energy, 



via the ansatz 

M*) = e^+^Mt) . (C6) 

Now one assumes that a weakly time dependent 
i/> a (t') can be replaced by the transition amplitude at 
the actual time t, ip a {t), thus neglecting the mem- 
ory effects of tp a . Using the integral relation (|C5|) 
again, one arrives at a expression similar to Eq. @ 

but with g (hu> - e h a — i" t + (e^J* - e^J replaced by 
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